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Moire patterns on Highly Oriented Pyrolytic Graphite surfaces due to dislocated graphene lay¬ 
ers were studied. We observed that the apparent corrugations of the moire patterns in scanning 
tunnelling microscopy images change as a function of tunnel junction bias. A simple geometric 
investigation of the atomic structure of the graphene layers generating moire patterns revealed that 
different atomic arrangements due to different twist angles can result in similar geometric moire 
periods such that, only a small fraction of the observed moire periodicities may coincide with the 
real atomic periodicity generating the moire system. Ab initio calculations showed that the band 
structure of moire patterns exhibit the fingerprints of those of bilayer graphene system preserving 
the Dirac cone. Moreover, our calculations in view of the correct atomistic modelling of the moire 
patterns showed that van Hove singularities in twisted bilayer graphene system with varying angles 
have different origins in their respective band structures. Consequently, our results shed light on the 
graphene like behaviour of the top most dislocated graphene layer on HOPG surfaces by showing 
that the layer does not act alone but the final graphene bilayer has electronic properties partially 
resembling graphene. 
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I. INTRODUCTION 

Highly Oriented Pyrolytic Graphite (HOPG) can be 
considered as the drosophila of surface science, espe¬ 
cially for the scanning probes. HOPG is a quite pop¬ 
ular crystal in experimental studies especially in sur¬ 
face science because of its smooth and chemically in¬ 
ert surface. It is widely employed as a calibration sam¬ 
ple for scanning tunnelling microscopy (STM) studies, 
as well as for educational purposes^^. HOPG is com¬ 
posed of Bernal stackecpHSl two dimensional hexagonal 
lattices, with carbon atoms at the lattice sites, popularly 
named as graphene layer^^M^. Nearest neighbour carbon 
atoms in a layer are bonded strongly with each other, 
while the l ayers ar e coupled only by weak van der Waals 
interact ion j ^ Q F U. 

Supe r pe riodic features were observed on HOPG 
surfacepn^HHI. Due to the weak van der Waals bond¬ 
ing between the layers of HOPG, topmost layer may be 
shifted or rotated by mechanical or chemical means. As 
a result of the rotation of the top layer, super-periodic 
structures called as moire patterns form on HOPGP®^ 
surfaces. Moire patterns occur in nature and they are 
fascinating phenomena already optically, and when it 
comes to such structures at the atomic leve l, nnder stand- 
ing their origin becomes even more tedioufl^EHUIl. Since 
the first observation of moire patterns on HOPG sur faces 
with STM, their origin is still a matter of debat^^^. 
Despite the fact that these formations were investigated 
in numerous studies, they are rediscovered in graphene 
researclP^M^. Due to its inevitable relevance to graphene 


based devices the importance of dislocated graphene on 
other gr aphen e layers on various systems drew serious 
attentiorP^^^^W^. Still the answers to the questions like 
’’When does the top most layer behave like graphene?”, 
or ’’Will that layer have other electronic properties” re¬ 
main yet to be settled. Localized effects also have to 
be understood and utilized or eliminated for any device 
application. 

By means of A ngle Resolved Photoemission Spec- 
trosc opy ( ARPES p^ * ^^ *, Landau Level Sp ectro scopy 
(LLS)P^tt 27 ] Quantum Hall Effect (QHE)P^^ mea¬ 
surements the graphene like behaviour of top most layers 
were studied. The electronic structure of the moire pat¬ 
terns and the van Hove (VHV) singularities observed on 
these systems were also i ntensively studied both theoret¬ 
ically and experiment ally^^^^ZElMTl ARPES results sug¬ 
gest that the top most layers of Graphene/SiC system al¬ 
ways behave like graphen^^ while the observed VHV sin¬ 
gularities and LLS measurements sugg est otherw ise, by 
simply showing a rotational dependenc d^^ * ^^ * ^ -^^. More¬ 
over, the QHE measurements support LLS and VHV 
measurements. Although moire patterns and their ge¬ 
ometric structure were considered in all studies, effect of 
the real atomic conformation of the moire patterns was 
not discussed in detail. The reason of the VHV singulari¬ 
ties from the band structures of the moire patterns is not 
complete yet, either. The lack of any moire structure ob¬ 
servation on HOPG by atomic force microscopy (AEM) 
is yet to be understood. So the question remains: ’’Are 
these structures comparable to the optical moire patterns 
or are they an outcome of the electronic structure of these 
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systems?” The variation of the apparent corrugations of 
the moire patterns on HOPG surfaces observed in STM 
measurements as a function of tunnel junction parame¬ 
ters were not investigated. In almost all the studies the 
corrugation was reported to be constant as a function of 
tunnel junction bia^^but for on^^, where the change of 
moire corrugation as a function of distance from defects, 
in resamblance to Friedel oscillations, was discussed. In- 
existance of the corrugation variation as a function of 
bias voltage was reported for moire patterns observed on 
Sid^ as well. 

Besides the experimental observations on HOPG and 
Graphene/HOPG systems, understanding these data be¬ 
came a challenging issue for theoretical studies. For in¬ 
stance, while the inspection of the HOPG surface be¬ 
came a benchmark for STM studieJ^, explaining the 
surprising triangular structures in spite of the hexago¬ 
nal lattice of graphite, which has two atoms in the sur¬ 
face unit cell, became a formidable challenge. Never¬ 
theless, this was successfully answered by realizing that 
the b site carbon atom which sits at the center of un¬ 
derlying hexagonal carbon ring of the second layer has 
slightly higher local density of states (LDOS) around 
Fermi energy compared to the carbon atom s itting at 
the a site (see Figure [^) of AB stacked graphit^^^^^EIl. 
Gonsequently, the first attempts to describe the reported 
STM images of moire superstructures were based on this 
facfP^, which had foreseen completely opposite bright re¬ 
gions. This was corrected by the model developed from 
the ab initio calculation for the DOS of AA, AB, and 
ABC stacked graphiteThis successful latter model 
was extended to investigate in detail the moire patterns 
on graphit^^. Obviously, all these approaches rely on 
perfect stacking results without performing a full Density 
Functional Theory (DFT) calculation on the rotated su¬ 
perstructures. Later on, some calculated STM images us¬ 
ing different DFT based methods were reported on moire 
superst ructure s due to twisted graphene layers on HOPG 
surfaceJS^^^^^. 

In this paper we report on the atomic structure of 
the moire patterns on HOPG surfaces and we unequiv¬ 
ocally show that many of the similar appearing moire 
patterns in STM images do not have the same atomic 
sub-structure. We report the apparent height corruga¬ 
tion variations of moire patterns as a function of tunnel 
junction bias measured with STM on HOPG surfaces for 
the first time. Our STS measurements and DOS calcu¬ 
lations show that moire structures behave more metallic 
than HOPG, but not graphene like. Our results show 
that the information on the real atomic structure lead¬ 
ing to each moire pattern is clearly invaluable for under¬ 
standing the electronic structure of this system, as they 
show for the first time that van Hove singularities ob¬ 
served on different moire systems do not have the same 
origin in their respective band structure. 


II. EXPERIMENTAL: PREPARATION AND 

MEASUREMENT OF MOIRE PATTERNS ON 
HOPG 

Various methods have been proposed for the repa¬ 
ration of moire patterns on HOPG surfaces^^^^ be¬ 
sides the random encounters on HOPG as well as on 
other graphitic systems like SiCP^ or multilayer graphene 
grown by Chemical Vapor Deposition (CVD) and trans¬ 
ferred on to other substrate^. Reported chemical meth¬ 
ods require time and extended processes, by means of 
which the HOPG crystal may be damaged permanently. 
We have prepared HOPG surfaces with moire patterns on 
them by means of a very simple and repeatable method. 
We have discovered that simply by drop casting cyclo¬ 
hexane on to HOPG surfaces (SPI^^ grade 2 and 3) and 
letting the samples dry in a clean environment, many dif¬ 
ferent moire structures can be obtainecP^. Pure and fresh 
(VWR^^ cyclohexane, HPLC grade) is used for each 
sample preparation. HOPG samples were thoroughly in¬ 
vestigated by STM prior to cyclohexane drop casting. 
After the drop cast and drying process, we have checked 
all the surfaces by means of micro Raman spectroscopy 
for probable contamination or residue from the cyclo¬ 
hexane drop casting and drying method. There were no 
observable remnants. Rarely (1 in 50 tunnelling zones) 
native moire patterns were observed on both pristine SPI 
grade-2 and grade-3 HOPG crystals. All the measure¬ 
ments reported here were performed under ambient con¬ 
ditions by use of a Nanosurf Easyscan2^^ STM and me¬ 
chanically cut PtsoIr 2 o tips. The experimental condi¬ 
tions were carefully controlled and a consistency check¬ 
list was developecP^. Many different moire samples were 
prepared and investigated by use of different tips. On 
the same sample reliable and consistent measurements 
for longer than half a year were possible. 

Since the site a atoms of the HOPG surface have 
less charge density compared to site b atoms (Fig- 
ure[^) STM gives triangular-atomic resolution on HOPG 
surfaced, where only the b atoms were observed. In some 
cases by use of quite high tunnelling current measure- 
ments, a atoms can also be observed by low temperature 
STIVP^. Intercalation of chemicals may result in hexago¬ 
nal resolutioiP^ on HOPG surfaces as well. On the moire 
patterned graphene/HOPG samples we always obtained 
triangular resolution on pristine HOPG terraces just next 
to the moire zones. More importantly, as shown in Fig¬ 
ure!^ and d, it was possible to achieve atomic resolution 
on moire patterns, simultaneously with the moire struc¬ 
ture observation. This was possible within a broad bias 
range (from 50 mV up to 800 meV). Still it was only 
possible to observe triangular atomic resolution on pris¬ 
tine HOPG terrace just next to the moire zone, only by 
50 meV, with the same tip during the same measure¬ 
ment. On clean monolayer graphene, transferred on to 
dielectric substrates, it was possible to o btain h exagonal- 
atomic resolution in STM measurement^^^^. The lack 
of total hexagonal resolution, but rather a partial hexag- 
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onal resolution on moire patterns indicate that top layer 
does not behave like graphene (Figure and d). Rather 
it must be the property of the top most layers resulting 
in the moire patterns. 


III. THEORETICAL 
A. Methods 

Efforts for understanding such intricate nature of the 
moire structures by electronic structure calculations of 
the twisted graphene layers for different rotation angles 
were performed using different methods. The major diffi¬ 
culty of computational studies was the number of atoms 
within the unit cell, which depends on the rotation an¬ 
gle and increases enormously with decreasing angl^^. 
Because of this, several studies were based on contin¬ 
uum theory, which exploited the effective Hamiltonian 
describing the l inear dispersion of Dirac cones of low en¬ 
ergy excitation j^^ l ^^^ ^. There were tight binding calcu¬ 
lations for the rotations, which can be describe d by a 
few thousands atoms within the unit cel P^ * ^^ * ^ ^. Accu¬ 
rate but computationally expensive ab initio calculations 
were very limi ted and only reported for relatively large 
rotation anglej^ ^^ l ^^ l Three different regimes emerged 
from these studie^^^^^ (1) For large rotation angles, two 
layers decouple from each other, so the linearly dispersed 
bands of Dirac cones of top and bottom layer remain 
unchangecP^^^. (2) In the intermediate regime, for ro¬ 
tation angles between 1° and 15°, Fermi velocity renor¬ 
malizes since the slope of linear Dirac bands decreases 
with decreasing rotation angle. Furthermore, the Dirac 
cone of the top layer and the rotated cone of the second 
layer intersect at finite energy so van Hove singularities 
emerge from these ” hybridiz ed” bands and clearly they 
depend on the rotation angl d^^ l ^^ l ^^ l ^^ l (3) Flat bands 
and vanishing velocities were reporte d for magical angles 
for rotations smaller than io| 3i | 37 | 5q 

We used VASP for the first-principle plane-wave cal¬ 
culations that employ density funct ional theory by the 
projected-augmented-wave potentiald^^^^. GGA is used 
to express the exchange-correlation potential and plane- 
wave cut-off energy is set to 500 eV. The Brilluoin zone 
integrations were carried out by using k-point meshes 
based on Monkhorst-Pack scheme defined according to 
size of the system (moire) unit cell, and the smearing 
parameter was set as 0.08 eV. The energy convergence 
was within 10“^ eV accuracy while relaxing the systems 
electronically. In view of our calculations here we discuss 
the properties of moire superstructures from the band 
structure and total density of states (DOS) as well as the 
STM images calculated by integrating the local density of 
states from the bias voltage (W) to th e Ferm i energy (Ef) 
as suggested by Tersoff and HamanrPSEIl^ Note that, it 
is necessary to use extremely dense k-point meshes in 
DOS calculations in order to avoid artificial oscillations 
and clearly show the features of van Hove singularities 



FIG. 1. (a) HOPG model; a type atoms are shown in blue 
and b atoms are indicated with red. (b) Atomic resolution 
STM image on HOPG, only b atoms are visible by STM {It= 
0.5 nA, Vb= 50 mV (sample bias)), (c) and (d) are STM 
images of different moire patterns. Mark the atomic resolu¬ 
tion on them, (c) It= 0.45 nA, Vb= -200 mV(sample), ob¬ 
served GMP=2.54 nm (d) /t=0.4 nA, V6=-50 mV (sample), 
GMP=1.6 nm. (e) Generation of the moire pattern (see text 
for details), (f) Generated moire pattern pi = l and gi = ll 
(Rotation angle=8.61°, GMP=1.64 nm and RMP=2.84 nm). 
Inset image shows calculated charge densities for this pat¬ 
tern. Observe the difference in the centers of super periodic 
formations. 


emerging in the DOS of the moire systems. In order to 
understand the effects of bottom layers on the charge 
density profile of the top layer we performed some pre¬ 
calculations on AB stacked graphite. We increased the 
number of layers from two to six gradually and observed 
no significant difference in the charge density profile of 
the top layer around the Fermi level. Therefore, we de¬ 
scribe moire patterns by twisting a bilayer graphene sys- 
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tern. We also performed tight binding (TB) calculations 
using the well-tested coupling parameters^!! in order to 
point out the critical differences with respect to ab initio 
calculations. 


B. Modelling the moire patterns: Real Moire 
Periodicity vs. Geometric Moire Periodicity 

We apply a practical scheme to describe the moire 
superstructure unit cell of twisted bilayer graphene as 
follows: Consider A A rather than AB stacked bilayer 
graphene as shown in Figure where only the top layer 
atoms of the model are visible. We ca n w rite any lattice 
translation vector of the top layer as 
by using the lattice vectors ^ and , where pi and 
gi are integers among themselves. In order to gener¬ 
ate the unit cell of twisted bilayer graphene, we need to 
find the corresponding lattice translation vector ^ of the 
bottom layer. ^ can be defined as the mirror symme¬ 
try of ^ with respect to p=0 line. Hence, ^ becomes 
—+ (gi +pi)^ . The angle between the p=0 line and 

^ and the angle between the p=0 line and ^ are same. 
Therefore, these two vectors coincide with each other on 
p=0 line by rotating the two layers towards this line by 0. 
(The angle of rotation of the top layer with respect to the 
lower layer within this scheme can be defined in terms of 
Pi and qi as 0rotation = 20 = 2arctan{\/3pi/2qi).) This 
rotated vector points to a lattice point of the twisted bi¬ 
layer. The other lattice vector of this two dimensional 
lattice of twisted bilayer graphene can be calculated by 
rotating this lattice vector by 7r/3 or 27r/3. These two 
vectors define the primitive unit cell of moire superstruc¬ 
ture unless gi-pi=3k, k being an integer. If that is the 
case, the primiti ve^ unit cell can be defined from a shorter 
lattice vector as R 2 = ((2pi+gi)^/3)-f((gi—pi)^/3) of 
which length is \Ri\/V3. Consequently, when pi is equal 
to one, 1/3 of all cases are special where real and observed 
(geometric) moire periodicities (RMPs and GMPs) are 
equal, while for the rest real moire periodicity is \/3 
times larger than the observed super-periodic structure 
(see Figure and Figure ©• If Pi is not equal to one, 
then there are even less cases where GMPs are the same 
as RMPs. We have performed first principles calculations 
for Pi = 1 and qi = 5, 7,10,11,12,16 cases. 

In Figure the rotation angle of the graphene layer 
is 8.61° (pi=l and gi=ll). Upon investigating the cen¬ 
tres of the round formations in such figures, they are ob¬ 
served to differ among themselves. Here the definitions 
of Real Moire Periodicity (RMP) and Geometric Moire 
Periodicity (GMP) become clear. In Figure]^ and b the 
calculated charge densities for two different systems, one 
being an RMP is different than GMP case and the other 
being a GMP is the same with RMP case, are given. In 
STM images such minute differences are almost impos¬ 
sible to identify and even in the model this small detail 
is easy to miss. In the literature moire periodicities are 



FIG. 2. (a) Calculated charge density for 0rotation=S.61°, 

GMP = 1.64 nm and RMP = 2.84 nm. Here D is the GMP 
but RMP is \/3D. (b) Charge density plot of a case where 
RMP=GMP=10.7 nm and Orotation =13.17°. In this case the 
moire unit cell is much smaller and consequently the systems 
has far less atoms to consider in ab initio calculations. 


assumed to follow the relation: D=a/(2sin(6>/2)) (where 
a is the lattice parameter) for bilayer graphene system. 
However, this relation only indicates the distance (D) be¬ 
tween the two bright spots observed, which is the GMP 
(see Figure Consequently, it would not be possible 
to deduct the real rotation angle using only the moire 
periods resolved in STM images. In such a case one can¬ 
not expect similar electronic structure for two dissimilar 
atomic configurations. Still the DOS of these systems 
show some surprising similarities, like a dependence of 
the energetic positions of van Hove (VHV) singularities 
on the rotation angle, which is a really lucky coincidence, 
which will be discussed below. 

We present two examples on the differences and sim¬ 
ilarities between RMPs and GMPs: When pi=l and 
gi=7 and pi=2 and gi=15 cases are considered, their 
GMPs are 10.72 Aand 11.4 Arespectively, which is pretty 
close and almost impossible to identify from STM im¬ 
ages. However, their RMPs are 10.72 Aand 39.59 
Arespectively, which shows the drastic difference between 
the amounts of atoms to be considered in the unit cells of 
both systems. When, the rotation angles for these cases 
are compared, the required angle for pi=l and qi=7 case 
is 13.17°, whereas this is 12.35° for pi=2 and gi=15. Due 
to such minute differences in angles, and since in an STM 
data only GMPs can be clearly identified, the real value 
of the rotation angle may not be evaluated from STM 
data even if atomic resolution was achieved. 


IV. RESULTS AND DISCUSSIONS 

A. Corrugation variation on moire patterns as a 
function of tunnel junction bias 

We collected STM images on various moire patterns 
with varying tunnel junction biases while the tunnelling 
current was kept constant. During the measurements 
moire structures were observed to be structurally stable. 
It was possible to observe atomic resolution on the moire 
patterns and investigate both atomic corrugation along 
with the moire corrugation by line scans as presented in 
Figure For each moire zone many STM images were 
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FIG. 3. 6.4 nm x 6.4 nm area scanned with 14= -45 mV 

(sample bias), It— 0.5 nA. Here the line scan on the moire 
pattern shows the atomic corrugation along with the moire 
corrugation 


obtained at different bias values while keeping the tun¬ 
nelling current constant. Than on every image several 
line scans were taken along 3 symmetry directions. From 
these data corrugation values were obtained. This was 
performed at least on ten different positions on the same 
image and on all the moire structures at every bias volt¬ 
age measured. Various tunnel junction biases were used 
from 50 mV to 800 mV without damaging the moire 
structures or the STM tip. In Figure lit we show the 
variation of the apparent moire corrugations observed on 
three different moire patterns (with different GMPs) as 
a function of the tunnel junction bias. In Figure |^-g we 
show that the atomic resolution was observable in a broad 
range of biases on this moire pattern. Decreasing corru¬ 
gation with increasing absolute value of tunnel junction 
bias was observed for these three super-periodic struc¬ 
tures. Measured GMPs of these structures were 2.54 nm, 
3.45 nm and 12 nm. While observed corrugation of two 
samples decreased with increasing tunnel junction bias, 
the third sample did not only show a decrease but also a 
varying behaviour in its apparent corrugation (both in¬ 
creasing and decreasing) as a function of increasing bias. 
In order to investigate the physical meaning of this ap¬ 
parent corrugation variation we have made a comparative 
charge density plot of the pi=l ^i=10 moire system with 
RMP=GMP=1.5 nm, as a function of different energy 
windows from Ef. These can be compared to different tip 
biases in STM measurements. Upon careful investigation 
of the local charge density distribution on the surface for 
the given energy window (i.e. Tersof-Haman calculated 
STM images), it is apparent that relative size of some of 
the charge density contours increase and than decrease 
with increasing bias (Figure [^,j and k). It is clear that 
STM tip cannot attain such a high lateral resolution but 
rather an averaged out image was obtained. Still, such 
variations in the electronic structure can be observed as 
corrugation variations in STM images. Both our exper¬ 
imental and theoretical data indicate that the apparent 
corrugation of the super periodic structures (i.e. moire 
patterns) varies with varying bias and this variance does 
not necessarily follow a decay. Such STM data and ab 
initio calculations show that the top most layer does not 
behave like graphene but the top two layers that generate 
the moire patterns on the HOPG surface have a new elec¬ 


tronic structure and the resulting surface may partially 
resemble graphene. 


B. Density of states of the moire patterns 

We performed scanning tunnelling spectroscopy (STS) 
measurements on moire patterns as well as on clean 
HOPG terraces just next to the moire patterned zones 
(Figure [^. The measured I-V curves were used for nu¬ 
merically calculating dln(I)/dln(V) vs. V curves which 
give the DOS of the surface of interest as commonly per¬ 
formed in the literatur^^. The immediate observation 
in the STS data is the relative metallic behaviour of the 
moire zones compared to bare HOPG terraces as seen 
in Figure and c. When the lineshapes of the STS 
curves are investigated, LDOS measurements of moire 
areas have a larger curvature (more parabolic than lin¬ 
ear) compared to the LDOS of HOPG terraces. This 
is an indicator of more metallic behaviour of the moire 
zones compared to the HOPG terraces. Another impor¬ 
tant observation is the sharp van Hove like states in the 
calculated DOS of the moire patterns, w hich were dis- 
cussed intensively in the recent literatnr d^^ l ^^ l ^^ l ^^ l ^^ l ^^ l 
. Several tight binding results were used to understand 
the nature of such features. According to previous re¬ 
ports energy values of the VHV singularities must follow 
a regular pattern as a function of the rotation angle of 
the top most layer. However, existence of a clear differ¬ 
ence between GMPs and RMPs suggest a more complex 
nature and the answer to this complexity lies in the band 
structure of each moire structure. 

In Figure and e we compare the band structures 
of the same moire pattern calculated using two different 
methods, ab initio planewave pseudopotential and tight 
binding calculations, respectively. These energy band 
structures are presented along the symmetry lines be¬ 
tween the high symmetry points of the Brillouin zone 
(BZ) of the moire superstructure unit cell of twisted bi¬ 
layer graphene as indicated in Figure [^. General fea¬ 
tures of moire band structures can be summarized as 
follows: There are linearly dispersed double degenerate 
bands around the K-point of moire BZ. However, the 
slope of these bands, which sets the Fermi velocity, de¬ 
creases with de creasi ng twist angle as consistent with the 
previous studie^^^^^. The interaction of the twisted top 
layer with its bottom layer does not open a band gap. 
Similarly, direct band gap at F point also depends on the 
twist angle. Interestingly, there are almost linear and 
fiat bands between the F and M points only in the ab 
initio calculated band structures. Similarly, some critical 
bands, as we will discuss below, between M and M' points 
as well as between K and F points are split in ab initio 
calculations while they are degenerate in tight binding 
results. Since the number of atoms within the moire su¬ 
perstructure unit cell depends on the twist angle and in¬ 
creases enormously with decreasing rotation angle, in the 
literature most of the studies were based on tight bind- 
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FIG. 4. (a) Bias dependent z-cormgation change of three different moire patterns, of GMPs 3.45 nm, 12 nm and 2.54 nm. 
Apparent z-corrugation of GMP=2.54 nm decreases with increasing bias while that of moire with 3.45 nm GMP does not 
exhibit a monotonic decrease and GMP=12 nm tends to stay almost constant in the same voltage range. Tunneling currents 
used in the measurements of 3.45 nm moire is 0.5 nA, 12 nm is 0.7 nA, and 2.54 nm moire is 0.5 nA. (b) During one of the 
scans on GMP=3.45 nm moire structure, bias was changed at several lines and moire corrugation was observed to decrease 
with increasing bias (It=0.5 nA, Vb was changed from -50 mV to -800 mV during scanning), (c) V6=-50 mV, (d) V6=-100 mV, 
(e) V6=-150 mV, (f) '\4=-200 mV , (g) V6=-250 mV (/t=0.5 nA for (c), (d), (e), (f), (g)). (h) Gharge density contour plots 
of pi = l gi = 10 at -400 meV. (i) zoom on to a small area - 400 meV (j) -800 meV (k) -1500 meV. Follow the charge density 
contours in the circled areas to observe the apparent size changes. 


ing calculations for the rotations, which can be described 
by a few thousands atoms within the unit. However, as 
we noted above there are some critical differences be¬ 
tween the ab initio and tight binding results. These may 
happen because of two reasons: (1) self consistent cal¬ 
culations are important for correct description of charge 
density distribution and (2) the inter-layer interaction, 
so the inter-layer hopping term needs to be adjusted for 
tight binding calculations. We must once more point out 
that with the rotation of the top layer the continuously 
changing observable in the super structure is the GMP 
not the RMP. As a result the band structure calculations 
must be investigated by keeping the atomistic structure 
of these patterns in mind. 


C. Reason of van Hove singularities in the DOS of 
moire patterns from their respective band structures 


We focused on one of the special cases in which RMP 
and GMP coincides and compared it to one of the more 
common cases where RMP is \/3 times the GMP. Dif¬ 
ferent types of critical points exist in the dispersion of 
two-dimensional structures like graphene. When the crit¬ 
ical point is a minimum or a maximum, a step function 
singularity occurs in the density of states, while at the 
saddle points, a logarithmic singularity arises. Hence, in 
general, saddle points are responsible for sharp peaks in 
the DO^. A nice example of this van Hove (VHV) sin¬ 
gularity is the sharp peaks around +1.5 eV and -2 eV 
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(a) 





(d) 



FIG. 5. (a) Calculated and measured DOS plots of HOPG. Size of the insets are 1.6 nm x 1.6 nm. 7^=0.5 uA and 14= -50 mV 
(sample bias) (b) Calculated DOS plots of HOPG, pi = l gi=ll moire and pi = l gi = 12 moire, (c) Measured DOS of two moire 
samples and HOPG. GMP of the first moire is 1.8 nm and the second moire is 1.6 nm. The size of the inset images are 12.8 
nm X 12.8 nm. Tunneling parameters; 7^=0.4 nA and 14=-50 mV . Mark that moire zones behave metallic compared to the 
HOPG surface, (d) Calculated band structure and density of states of pi = l qi—11 moire structure with ab initio planewave 
pseudo potential method, (e) Calculated band structure of pi=l gi = ll moire structure with tight binding method, (f) High 
symmetry points used for the calculation of band structures inside the moire Brillouin zone. High symmetry points are marked 
as a guide to the eye. 


at DOS of graphene (or graphite), which are due to the 
saddle points at M point of its BZ. In Figurethe band 
structure of pi=l and gi=10 moire is shown, which is 
one of the special cases where RMP is the same as GMP. 
In this system the rotation of the layers with respect to 
each other is 9.43° and the moire period is 1.496 nm. 
When we inspect the band structure and corresponding 
DOS thoroughly, it is evident that the VHV originates 
from the saddle point of band edges between F and M 
points, and in this case it is closer to the M point. We 
also show and discuss this in Figure In Figure 

we present the band structure for one of the common 
cases, namely pi=l and ^^ 1 = 12 , which has a rotation an¬ 
gle of 7.93° with GMP of 1.78 nm and RMP of 3.082 nm. 
Again upon inspection of the band structure we see that 
the VHV of this case originate from the saddle point of 
the bands along the line connecting F and K points. We 
can mark this between the M and M' points as well (see 
Figureand f). This is better visible in Figure [^. 

From our results it is clear that the VHV singularities 
occurring in the DOS of different moire patterns do not 


have the same origin in their band structure and this is 
clearly related to the existence of GMPs and RMPs. So 
far this has not been recognized in the reported stud¬ 
ies. Rather, the emergences of VHV singularities were 
associated to intersection of Dirac cones of top and bot¬ 
tom layer exactly at M point of moire BZ. In Figure |^, 
we present the iso-surface of the 7r-band of the graphene 
band structure in the irreducible part of the graphene 
BZ around the K point. Essentially, this is a cosine band 
within tight binding approach and the linear part around 
the K point is known as Dirac cone. However, the Dirac 
cone diverges from a circular cone as the energy increases. 
We want to emphasize that the energies at which the 
VHV singularities occur are not close to the Fermi en¬ 
ergy. On this band diagram, we highlighted the energies 
corresponding to the VHV singularities of the twisted 
bilayer graphene. The inner most curve corresponds to 
100 meV. At this low energy the cross section is circular, 
which is in the linear part of the Dirac cone. Apart from 
these cases this constant energy cross section of 7r-band 
of graphene is more triangular, which becomes clear at 
























FIG. 6. Calculated band structure (ab initio planewave pseudopotential method) and DOS of (a) pi = l and ^i=10 moire 
structure, (9.43°, special case of 1/3, where RMP=GMP=14.96 A) has VHV at about 0.58 eV. (b) pi=l and gi = 12 moire 
structure, (7.93°, GMP=17.8 Aand RMP=30.82 A) has VHV at about 0.46 eV. (c) Energy band diagram of tt band of graphene 
within the irreducible BZ. Energies corresponding to the VHV singularities of various moire structures are indicated, (d) The 
relation of bilayer graphene and moire BZs at the emergence of Van Hove singularity of pi=l gi = 10, and (e) of pi = l gi = ll 
moire structures, (see the supplementary video for further information.) 


higher energies. Apparently, these are outside the low 
excitation region where the continuum theories are ap¬ 
plicable. If we had linear Dirac cones, then the inter¬ 
section of these cones from top and bottom layer should 
occur exactly at M point of moire BZ. As more of such 
non-circular cones intersect each other, the intersection 
generating the VHV singularities diverges form linearity 
(M point is a special case). The cone intersection does 
not occur as commonly represented for bilayer graphene 
at small energy values around the K point. Consequently 
the cones emanating from the K points do not intersect 
at the same places in the BZ for every rotation angle, 
as exemplified in Figure and e. Furthermore, this in¬ 
tersection point depends on the orientation of the moire 
BZ, so the origin of VHV singularity in the special cases 
where GMP=RMP and general cases where they are not 
equal are completely different. 


V. CONCLUSIONS 


Our combined STM and ab initio studies on the moire 
patterns on HOPG surfaces clearly showed the impor¬ 
tance of proper atomistic modelling of the moire patterns 
on such systems. Although DOS of RMPs and GMPs 
might appear similar, the band structures of these sys¬ 
tems differed in detail. The emerging VHV singularities 
in their respective DOS had completely different origins 
from the energy bands of corresponding moire systems. 
We have shown that the corrugations of different moire 
patterns change as a function of the tunnel junction bias 
differently from each other in STM measurements. Al¬ 
though STM may not have the necessary spatial resolu¬ 
tion to image the localized charge distribution of these 
moire structures, we have evidence, both theoretically 
and experimentally, that the corrugation of the moire 
patterns do not change monotonically. While we show 
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that the corrugation of the moire structures is not con¬ 
stant as a function of energy, we also show evidence for 
the non-graphene like behaviour of the top most layer of 
the moire structures. Rather, we have shown that the 
graphene bilayer system generating the moire structures 
behave as a large molecule. 

In this study we have shown a very simple, reliable 
and repeatable method of producing moire patterns on 
HOPG surfaces without damaging the crystals. We have 
clear evidence that the moire/HOPG systems are perfect 
for being studied under ambient conditions with STM. 
Finally our results indicate that although moire patterns 
on HOPG surfaces have mainly an electronic origin, the 
corrugation measurements show the hint of a minute 
structural origin, which is yet to be tested. 

Further measurements like dl/dV maps of the moire 
patterns with high corrugation variance will also help un¬ 
derstand the localized electronic structure of the moire 
patterns. The corrugation variation of moire systems are 
yet to be tested on SiG or graphene bilayers grown by 
GVD. Local electronic measurements on different moires 
is another challenge for device applications. Our results 


have also direct relevance to moire patterns observed on 
other layered systems such as hBN, or Graphene/hBN 
systems, especially in their modelling, which attracted 
recent attention^ . Gorrugation variation in such sys¬ 
tems is another issue that awaits testing. The ease in the 
preparation of the samples we used for this study may 
help many others to pursue similar methodologies in sim¬ 
ilar layered systems. Our data has also direct relevance 
to graphene device literature where the identification of 
mono/bi-layer or multilayer graphene systems is invalu¬ 
able while producing graphene based devices, which ide¬ 
ally function at room temperature and especially under 
ambient conditions, for real life applications. 
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